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Abstract 

Assuming a Friedmann universe which evolves with a power-law 
scale factor, a = t n , we analyse the phase space of the system of 
equations that describes a time- varying fine structure 'constant', a, 
in the Bekenstein-Sandvik-Barrow-Magueijo generalisation of general 
relativity. We have classified all the possible behaviours of ct(t) in 
ever-expanding universes with different n and find new exact solu- 
tions for ot(t). We find the attractors points in the phase space for 
all n. In general, a will be a non-decreasing function of time that 
increases logarithmically in time during a period when the expansion 
is dust dominated (n = 2/3), but becomes constant when n > 2/3. 
This includes the case of negative-curvature domination (n = 1). a 
also tends rapidly to a constant when the expansion scale factor in- 
creases exponentially. A general set of conditions is established for 
a to become asymptotically constant at late times in an expanding 
universe. 



1 Introduction 

Stimulated by the observations for small variations in atomic structure con- 
trolled by the fine structure constant in quasar absorption lines at redshifts 
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z = 1 — 3, []T|, 0, [J, there has been much recent interest in the theoretical 
predictions of gravity theories which extend general relativity to incorpo- 
rate space-time variations of the fine structure 'constant'. These have been 
primary formulated as Lagrangian theories with explicit variation of the ve- 
locity of light, c, f|, H, U, or of the charge on the electron, e, jlO], |, H, 0. 



Theories of the latter sort offer the possibility of matching the magnitude 
and trend of the quasar observations and have been studied numerically and 
by means of matched asymptotic approximations. They are also of partic- 
ular interest because they predict that violations of the weak equivalence 
principle should be observed at a level that is within about an order of mag- 
nitude of existing experimental bounds [|ll, 12, 13fl . They are consistent with 



all other astrophysical and experimental limits of time variation of the fine 
structure constant and predict effects of the microwave background radia- 
tion, primordial nucleosynthesis, and the Oklo natural reactor that are too 
small to conflict with current observational bounds. A range of variant the- 
ories have been investigated with attention to the possible particle physics 
motivations and consequences for systems of grand and partial unification in 



references |14| , [15], [L6[ . In this paper we will give a full qualitative analysis of 
the properties of Friedmann cosmological models in a sub-class of these the- 
ories developed initially by Bekenstein |17| to generalise Maxwell's equations 



to include varying e and then generalised by Sandvik, Barrow and Magueijo 
HlO| to include gravitation. We refer to these as BSBM theories. We provide 
a phase-space analysis of the non-linear propagation equation for the scalar 
field which carries the variations of the fine structure constant. Some new 
exact solutions are also given and all the asymptotic behaviours classified. 

2 The BSBM Theory 

2.1 The Cosmological Evolution Equations 

We will assume that the total action of the Universe is given by BSBM: 

S J d X\J Q (y£--grav ^matter ^-"t/j £ern& ^) > (1) 

In the BSBM varying a theory, the quantities c and h are taken to be con- 
stant, while e varies as a function of a real scalar field ip, with 
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e = eoe^ (2) 

where = ^d^d^i/;, u is a coupling constant, and C em = —\f^ u f liu . The 
gravitational Lagrangian is the usual C g = — j^^R, with R the curvature 
scalar and we have defined an auxiliary gauge potential a M = eA^, and field 
tensor f^ u = eF^ u = d^a u — <9„a M , so the covariant derivative takes the usual 
form, = + ieoa^. The dependence on e in the Lagrangian then occurs 
only in the kinetic term for e and in the F 2 = f 2 /e 2 term. 

The universe will be described by a homogeneous and isotropic Friedmann 
metric with expansion scale factor a(t) and curvature parameter k. Varying 
the total Lagrangian we obtain the Friedmann equation (G — c = 1) for a 
universe containing pressure-free matter and radiation: 

2 = Y ( Pm (1 + ICI e" 2 ^) + p r e- 2 * + p, + PA ) - A ( 3 ) 

where the cosmological vacuum energy p A is a constant given by A/ (87r), and 

/°V = f V* an d C — ^em/ p m is the fraction of the matter which carries electric 
or magnetic charge. 

For the scalar field we obtain the evolution equation 

^ + 3#V> = --exp[-2V>]Cp m (4) 

UJ 

where H = a/a is the Hubble parameter. The conservation equations for the 
non-interacting radiation and matter densities, p r and p m respectively, are: 

p m + 3Hp m = (5) 
p r + AHp r = 2^p r . (6) 

so p m oc a~ 3 . However, the last relation can be written as p r + 4Hp r = 0, 
with p r = p r e~ 2 ^ oc a -4 . 

The Friedmann models with varying a have been shown |J to have the 
property that when ( < the homogeneous motion of the ip does not create 
significant metric perturbations at late times. This means that far from 
the singularity we can safely assume that the expansion scale factor follows 
the form for the Friedmann universe containing the same fluid when a does 
not vary (£ = = ip). The behaviour of ip then follows from a solution 
of equation (^) in which a(t) has the form for a Friedmann universe for 
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matter with the same equation of state in general relativity when £ = = ip. 
This behaviour is natural. We would not expect that very small variations 
in the coupling to electromagnetically interacting matter would have large 
gravitational effects upon the expansion of the universe. Thus, in this paper 
we will provide a complete analysis of the behaviour of the solutions of the 
non-linear propagation equation (|) for appropriate behaviours of a(t). 



2.2 An Approximation Method 

We consider spatially flat universes (k = 0) and assume that the expansion 
scale factor is that of the Friedmann model containing a perfect fluid: 

a = t n (7) 

where n is a constant. The late stages of an open universe containing fluid 
with density p and pressure p obeying p+3p > can be studied by considering 
the case n — 1. We rewrite the wave equation (|4]) as 

(^a 3 ) = p m a 3 exp[-2^] 

Therefore, since p m a 3 is constant this reduces to a Liouville equation of the 
form 

(^a 3 )- = iVexp[-2V>] (8) 
where N is a constant, defined by 

N = -^ Pm a 3 . (9) 
to 

We shall consider first the cosmological models that arise when the defining 
constant N is negative. This arises when the constant ( < 0, indicating that 
the matter content of the universe is dominated by magnetic rather than 
electrostatic energy. The value of ( for baryonic and dark matter has been 
disputed |L2], [16, 10| . It is the difference between the percentage of mass in 



electrostatic and magnetostatic forms. As explained in JHJ, we can at most 
estimate this quantity for neutrons and protons, with ( n w ( p ~ 10~ 4 . We 
may expect that for baryonic matter ( ~ 10~ 4 , with composition-dependent 
variations of the same order. The value of ( for the dark matter, for all we 
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know, could be anything between -1 and 1. Superconducting cosmic strings, 
or magnetic monopoles, display a negative (, unlike more conventional dark 
matter. It is clear that the only way to obtain a cosmologically increasing 
a in BSBM is with £ < 0, i.e with unusual dark matter, in which magnetic 



energy dominates over electrostatic energy In []10[] we showed that fitting 
the Webb et al results requires Q/uj = —2 ± 1 x 10 -4 , where ( is weighted by 
the necessary fractions of dark and baryonic matter required by observations 
of the gravitational effects of dark matter and the calculations of Big Bang 
nucleosynthesis. We note also that in practice ( might display a significant 
spatial variation because of the change in the nature of the dominant form 
of dark matter over different length scales. For example, a magnetically 
dominated form of dark matter might contribute a negative value of ( on 
large scales while domination of the matter content by baryons on small 
scales would lead to ( > locally. We will not discuss the effects of such 
variations in this paper. 

2.3 The Validity of the Approximation 

We have assumed that the scale factor is given by the FRW model and 
then solved the ip evolution equation. This is a good approximation up to 
logarithmic corrections. Here is what happens to higher order. 
We take the leading order behaviour in (|3|) 

• \ 2 o 



, ,P m (l + IC|e- 2 ^) (10) 



Now if we take a = t 2 ^ 3 so 
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and solve equation (|J) we get asymptotically, 

if) = ^ln[2iVln£] -> \n[\nt] (11) 

to leading order at late times. Suppose we now re-solve (flDI ) with the \( \ e~ 2 ^ 
correction included 

°V K _i (1 + i C i^) K _i( 1 + ^_) (12) 
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Note that the kinetic term which we neglected is of order 

/ 2 1 
uoip ~ 



tfln'« <13) 
and so is smaller than the term we have retained. Solving equation (|II]) we 
have 



In a 



2|C| 



, , ,lnt /hit / lnt 



(14) 



Note that when |d — > this gives the usual In a = | Int. When |d 7^ we 
have 



2 (CI L r /lnt, lnt 



and 



t 2 / 3 (lnt) 



CI/3 



(15) 



where \(\ is small and so the corrections to the a = t 2//3 ansatz are small. In 
terms of the Hubble rate: 



3t 3tlnt 

If we include the kinetic corrections to equation ( pT0|) then as if) - 1 



2tlnt 



a / 9t 2 



(l + K|e- 2 ^) + — ^ «_(l + JiL + TX 



lnt \n 2 t 



(16) 



where 
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3ttlu 



So, if x = lnt, we have 
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In a 



y/x 2 + \C\x + S + S J 



dx 



x^/x 2 + |C| x + S 



2 



dx 



^x 2 + \(\x + S 



(17) 



■In a = \/x 2 + |C| x + S - In 



2S + |C| x + 2y/Sy/x 2 + \t\x + S 



2 



2 v / z 2 + \(\x + S + 2x + |C| 



;is) 



Again, as t — > oo the leading order behaviour is that found in equation (15). 
In the radiation era we have an exact solution of equation (B) with 



2i\) = ln(8N) + ilnt 
so the corrections to the Friedmann equation look like 

2 



At 2 



(1 + Kle- 2 ^) 



1 
it 2 " 



1 + 



ICI 



IN t 1 ' 2 



(19) 



(20) 



and these corrections fall off much faster than in the dust case. Again, our 
basic approximation method holds good to high accuracy. 



2.4 A Linearisation Instability 

Despite the robustness of the basic test-motion approximation that we are 
employing to analyse the evolution of ip(t) as the universe expands, there is a 
subtle feature the non-linear evolution equation ([8]) which must be noted in 
order that spurious conclusions are not drawn from an approximate analysis. 
We see that the right-hand side of equation ([8]) is always positive. Therefore 
ip can never experience an expansion maximum (where ip = and ip < 0) and 
therefore ip(t) can never oscillate. However, if we were to linearise equation 
(H), obtaining 

(^a 3 )- = iVexp[-2^] w N(l -2i/) + 0(tp 2 )) 
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then for ip > 1/2 the right-hand side takes negative values and pseudo- 
oscillatory solutions for ip would appear that are not the linearised approxi- 
mation to any true solution of the non-linear equation (|8|). Care must there- 
fore be taken to ensure that analytic approximations are not extended to 
large ip and that numerical analyses are not creating spurious spirals in the 
phase plane by virtue of a linearisation procedure; for a fuller discussion see 
ref. 0. 

These considerations can be taken further. It is possible for ip(t) to de- 
crease, reach a minimum and then increase. But it is not possible for ijj(t) 
to decrease if it has ever increased. A second interesting consequence of this 
feature of equation (||]) is that it holds true even if a(t) reaches an expansion 
maximum and begins to contract. Thus in a closed universe we expect ip and 
a to continue to increase slowly even after the universe begins to contract. 
This will have important consequences for the expected variation of ip and a 
in realistically inhomogeneous universes. 



3 Phase-plane Analysis 

3.1 A transformation of variables 

In this section we will look at the i/j equation of motion (|8])when the expansion 
scale factor takes the power-law form (|7|). The evolution equation for the field 
then becomes: 



d 
dt 

with N > 0. We introduce the following variables : 



(^jt 3n ^j = iVexp[-2?/>] (21) 



3n 1 

x = ln(t) y = ip- Ax- B A — I — B = - ln(iV) (22) 

and rewrite (H) as: 

y" + (3n-l)(y' + l-^)=e-*y, (23) 

where ' = d/dx. This second-order differential equation can be transformed 
into an autonomous system by defining v — ^- and u = y: 
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% = e-^ + (l-3n)(, + l-^) (24) 

du 
dx 

We see that the system has a finite critical point at 

1 Sr) 

(« e> « c ) = (--lB[(l-3n)(y-l)],0) (25) 

when n G (|, |) and it has a infinite critical point at (u c ,v c ) = (+oo,0) 
when n = | or n = |. The finite critical points correspond to the family 
of exact solutions of the form ip = B + A\n(t) found in ||. In the original 
variables these solutions are 

, L r 2N , /2-3n\ , , , 

ip = - In + In (t ). 

y 2 L (2-3n)(3n-l) J \ 2 J w 

In order to analyse the system fully, we will study finite critical points 
and the critical points at infinity separately. We will also distinguish several 
domains of behaviour for n: n = ^, n = |, n G (|, |), n < | and n > |. In 
some cases we will supplement our investigations by numerical study of the 
system 



3.2 The Finite Critical Points 
3.2.1 The Cases n G (|, |) 

In the domain where n G (3, |) the system ( p4[) does not appear to be exactly 
integratable. So in this section we will study the behaviour of the system 
near the critical points and at late times. 

There is a finite critical point (u c ,v c ) = (— |ln [(1 — 3n)(^ — 1)],0) for 
n G (3, |). Linearising the system ^ about it we obtain: 



dV 
dx 

— = V 

dx 



(3n-l)(3n-2)C^ + (l-3n)V r (26) 
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where V = v — v c and U = u — u c . Since the characteristic matrix is non 
singular the critical point is simple. Hence the non-linearised version of this 
system has the same phase portrait at the neighbourhood of the critical point. 

The eigenvalues ^ 2 and corresponding eigenvectors x^ 1,2 °f t ne system 
(|26|) are: 



1 - 3n ± V9 - 42n + 45n 2 , , 

£1,2 = o ( 27 ) 



, , 1 - 3n ± V9 - 42n + 45n 2 . 
X fl ' 2 = ( 2 ' 1} 

These eigenvalues will be complex numbers for n G (|, |) and they will be 
pure real numbers when n G (|, |). Notice that for both cases the real part of 
the eigenvalues is always negative, so the critical point is a stable attractor. 
The general solution of the linearised system (|26|) can be expressed as: 



(U(x),V(x)) = Ae^xi + Be^ x x 2 

where A, B are arbitrary constants and x = \n(t). From the transformations 
(H) we can obtain explicitly the expression for if)(t) near the critical point. 
There are three possible behaviours of the solutions near the critical point: 

3.2.2 Pseudo-Oscillatory ip behaviour 

When n G (|, |) the linearised evolution for the ip field exhibits damped 
oscillations about its asymptotic solution as t — > oo: 



where 



^ = i ln ' (l-3 n H3, i -2) l^< 2 - 3 """< t) < 28) 
+t a (5 cos[/3 ln(t)] + A sm[(3 \n(t)]) 



1 — 3n , , 

a = ^— (29) 

/J = V9 + 42 "- 45 " 2 (30) 
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Figure 1: Numerical plots of the phase space in the (w,v) coordinates, and the tp 
evolution with x = log(t) for n = ^ . The '+ ' sign is a saddle point and the triangle 
is a stable node. 

We note that this family of solutions and asymptotes includes the case of the 
radiation-dominated universe (n = 1/2). The phase portrait and the ip vs. 
hit for this case are shown in Figure 1. 

But in the case of a universe containing the balance of matter and ra- 
diation, displayed by our own, it need not be the case that the asymptotic 
behaviour, displayed by the exact solution for the critical point, is reached 
before the radiation-dominated expansion is replaced by matter-domination, 



see references JT0[ and [[| for further discussion of this point. However as we 
discussed above these oscillations are an artefact of the linearisation process 
and the part of the solution ( ^8|) controlled by the constants A and B is only 
valid for small times, hence we called this behaviour pseudo-oscillatory. 

3.2.3 Non-oscillatory behaviour 

When n £ (|, |), the ip field will approach its asymptotic behaviour in a 
non-oscillatory fashion as t — > oo: 



ib(t) = -ln[- ^ -] (31) 

^ W 2 L (l -3n)(3n-2) J V ' 

+ ^ (2 - 3n) ln(t) + t 5 ' 1 (A + Bt 2 ^) 
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where 



5 = 



7 



1 -3ra 



-V9-42n + 45n 2 



(32) 
(33) 



Note that 5 + 7 < in the domain n G (|, |), so at late times (t — > 00), 
both the pseudo-oscillatory and non-oscillatory solutions case will approach 
the asymptotic solution defined by the appropriate value of n. 

3.2.4 Intermediate behaviour 



A transition between the pseudo-oscillatory and non-oscillatory regimes hap- 
pens when (3 — and n — |; th 
the vicinity of the critical point: 



pens when (3 = and n = |; the ip field then has the following solution in 



^(t) = Un(N) + ±\n(t)+At- 



(34) 



The phase plane structure and evolution of ip vs. In t for this case are shown 
in Figure 2. 




Figure 2: Numerical plots of the phase space in the (w,v) coordinates, and the 
tp evolution with x = log(t) for n = |. The '+' sign is a saddle point and the 
triangle is a stable node. 
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3.2.5 Overview 



The late-time solution of the ip field solution is given by asymptotic behaviour 



of equations (|28Q or which is: 



V> -> - Inb T^r rl + - (2 - 3n) ln(t) 

^ 2 L (l-3n)(3n-2) J 2 V 7 w 

This shows that the solutions fl28|) and ([H]) generalise the ones found in ||, 
since they can be obtained setting A = B = in (|2S|). In particular, the case 
of a radiation-dominated Universe (n = |), we have from (f28|): 

= ~(log(8)+log(iV)+log(*)) (35) 

, ,_i / R r 5v / 31og(t) 1 , , . r 5yglog(t ) 1 

+t * i^cos +^4 sin 

14 4 

A full mathematical summary of the change in structure of the phase space 
with changing n for the system (||) is given in the Appendix. There we 
include cosmologically unphysical values of n and show how the critical point 
structure bifurcates with the change in value of n. 

3.3 The Critical Points at Infinity 

In order to describe the qualitative evolution of the system, we must de- 
termine the behaviour of the system (|24| ) near the critical point (u c , v c ) = 
(+oo, 0). In order to bring the critical point to a finite value, we define: 

u = — \n(w) 

Using the new coordinate w we can re-write the system ( |24"D as: 

dw , . 

— = -2wv (36) 
ax 

This system has critical points on the (w, v) plane when (0, —1 + ^y) and 
((1 — 3n) f ^ — 1J , 0). Note that the second critical point is just the same as 
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the one we have analysed in the previous section. In this subsection we will 
then only analyse the critical point (w c , v c ) = (0, — l + since it corresponds 
to the case where u — > +00. 

Proceeding as before, and linearising fl3~6|) about (w c ,v c ) = (0,-1 + ^) 
we obtain: 

(2 - 3n)W (37) 
W + (1 -3n)V 

where V = v — v c and W = w — w c . Again, the characteristic matrix of the 
system is non singular, and the critical point is simple. Hence, system Q3"7| ) 
will have the same phase portrait as (|36| ) in the neighbourhood of the critical 
point. 

The eigenvalues ^ 2 and corresponding eigenvectors x^ 1 ' 2 °f t ne system 
are: 

Ci = l-3n, X 5l = (0,l) (38) 

e 2 = 2-3n, X C2 = (1,1) (39) 

These eigenvalues are always real. The critical point is an attractive node 
for n > |, a saddle point when n G (|, |), and it will be an unstable point 
when n < |. The general solution of the system in the neighbourhood 
of the critical point (w c , u c ) = (0, —1 + ^), can be expressed as: 

(w(x) - w c , v(x) - v c ) = Ae^ x Xi + Be^ x x 2 

where A, B are arbitrary constants and x = ln(t). Therefore, near the critical 
point (u c , v c ) = (+oo,0): 

m = \H%) ( 4 o) 

so the scalar field is constant. 

Note that in the domain n G (|, |) this is just a transitory solution since 
the critical point is a saddle point. In the domain n < | it is an unstable 
critical point, possibly relevant as an early-time solution to braneworld cos- 
mologies in the high-density regime where the Hubble expansion rate of the 



dW 

dx 
dV 

dx 
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universe is linearly proportional to the density, so n = 1/4 for a braneworld 
containing radiation, n = 1/12 for a massless scalar field, and n = 1/6 for 

dust. However, in the t — > limit the assumption that the t/> and ^exp[— 2-0] 
terms can be neglected in the Friedmann equation (|3]) will break down. 

3.3.1 The cases of n > 2/3 and de Sitter expansion 

An interesting case exists when n > |, since the critical point is a stable at- 
tractor and so this means that the constant-'?/' behaviour (KH) is the late-time 
attractor, in agreement with the conclusions of references |10[ and M. This 
is an important feature of a universe which exhibits accelerated expansion in 
its late stages (n > 1). It means that the present value of a is the asymptotic 
one. It also means that variations in a are turned off by the domination of 
the expansion dynamics by negative curvature or by any quintessence field 



iOfl , ||. This property may provide important clues to explaining why our 
universe possesses small but finite curvature or quintessence energy today: if 
it did not then the fine structure constant would continue to increase until 
it was impossible for atoms and molecules to exist f|. 

In the n > 2/3 case we can find a detailed asymptotic solution for equation 



(|2T| ) which has the form 

v> = c + B( t + to r^ + Ne ^}~ 2 n c \ t + t ) 2 - 3n 

with C, to constants. We see immediately that for this solution ip approaches 
a constant as t — > oo. Of particular interest is the case of a curvature- 
dominated open universe, which has n = 1. The phase plane structure and 
the evolution of ip vs. Int for this case are shown in Figure 3. 

We see that this approach to constant behaviour occurs for all universes 
that accelerate (n > 1) and so we would expect to find it also in the case of 
a de Sitter background universe with 

a(t) = exp[At] 

where A > is constant. 

Substituting this in equation (§) we find a late-time asymptotic solution 

if> = C + Bexp[-3\t] - N ( 3Xt + X ) exp [-2C - 3At] -> C 

9A 
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Figure 3: Numerical plots of the phase space in the (w,v) coordinates, and the tp 
evolution with x = log{t) for n = 1. The '+ ' sign is a saddle point and the triangle 
is a stable node. 



as t — > oo. 

Notice that we were not able to fully describe the nature of this critical 
point in the cases where n = | or n = | since one of the eigenvalues of the 
systems is zero. In order to do so, we will study these two cases individually, 
in particular the case n — | is important since it describes the scale factor 
evolution on a dust dominated universe. 

3.3.2 The n = |, a{t) = U case: an exact solution 

The n = | case can be exactly integrated. It is of interest as an exact solution 
in its own right but it corresponds to the case of a universe whose expansion 
dynamics are dominated by the effects of a fluid with a p = p equation of 
state, or a massless scalar field. It also describes the behaviour of the ip field 
in an anisotropic universe of the simple Kasner type. We see when n = 1/3 
the system (|23j) has the form: 



^ = e- 2 " (41) 
ax 



du 
dx 

This integrates to give 
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v 2 + e -2u = E 2 

where E is a constant. Hence, we have two possible solution branches: v (u) = 
±V E — e~ 2u . Both, positive (+) and negative (— ) solutions for v lead to the 
same result. Choosing the positive branch of the v (u) solution and using 
equations (pLlD we obtain: 

u{x) = - \n(E) + lnfcosh (E (d - x))} 
where G\ is an integration constant. From (|22|) we have a solution for ip(t): 

*P(t) = -y^- - HE) + 2 ln (^) + He*P[ECi]r E + t E expf-^d]] 
and the asymptotic limit when t — > oo gives 

m^(\ + E) \n(t + t ). 

with to constant. The phase space and ip(t) vs. Int evolution is shown in 
Figure 4. 




x 



Figure 4: Numerical plots of the phase space in the (w, v) coordinates and time 
evolution of ip(x) n = |. The central line which passes through the origin is an 
attractor for w > 0, and an unstable line for w < 0. 
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3.3.3 The n = § , a = ts case 

The case of a dust-dominated universe is mathematically special because of 
the presence of a zero eigenvalue in the stability analysis. It is illuminating 
to consider this case separately with an asymptotic analysis that extends the 
earlier study in ||. Consider equation ( 2lj) with n = 2/3. If we introduce 
the new variable x = ln(£), then 

if/' + ij/ = N exp[-2ij}] (42) 
At large x, the asymptotic form of this equation has the form: 

V = g ln[27V(x + *„)] - ^ E ( JVx o) " + (43) 

where C and Xo > are arbitrary constants. The leading order behaviour as 
t — > oo is therefore (cf.flTTD) 

ip ~ -ln[2JVln(t)] 

Stability analysis of the n = 2/3 asymptote Performing the coordinate 
transformation it = — | ln(w) on the system (24) when n = 2/3, we have: 



— = w - v (44) 
ax 

— — = —2vw 

ax 

This system has a critical point at the origin of the (w,v) plane. It cor- 
responds to the case where u — ► oo. The linearisation about this critical 
point gives two eigenvalues: A x = —1 with the eigenvector \i — (0, 1) and 
A2 = with the eigenvector x 2 = (1,1) . Since we have a zero eigenvalue, the 
stability is determined by the non-linear behaviour and is one of Lyapunov's 



'critical' cases p7|,p8| |29|| . We apply a linear transformation to split the 
system into critical and non-critical variables; where the critical variables are 
those eigenvectors with zero eigenvalue and the non-critical variable are the 
others. Then we apply a non-linear transformation which will eliminate the 
influence of the critical variables upon the non-critical ones at the leading 
order. This gives (w,v) — > (W, W + V), and the system ( 44|) becomes: 
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W' = -2W(W + V) 
V' = -V + 2W(W + V) 



(45) 



with the critical point at (W,V) = (0,0). The Lyapunov procedure for the 
system (45) is to set the linearly stable variable, V, equal to zero in the W 
equation so, 



W' = -2W 2 



and so the second-order analysis shows that critical point (W, V) = (0, 0) is 
unstable. In fact this critical point is a saddle point as can be seen from 
the numerical-phase-plane plot figure [5]. The unstable part correspond to a 
non-physical range of the cosmological variables, since it gives w < 0. The 
stable part corresponds to the range of values where w > 0. With respect 
to our (approximate) exact solution, the unphysical case will correspond to 
the range of x < 0, since they lead to a < 0. Hence the asymptotic solution 
fl43| ) is the stable late-time behaviour of the dust universes with small (. The 
phase plane structure and the evolution of ip vs. lni is shown in Figure 5. 




Figure 5: Numerical plots of the phase space in the (w,v) coordinates, and the ip 
evolution with x = log{t) for n = |. The triangle is a saddle point for w < and 
a stable node for w > 0. 
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4 Some general asymptotic features 



4.1 Models with asymptotically constant ip and a 

We have seen that ip, and hence the fine structure constant, a, tends to a 
constant at late times in accelerating universes with power-law and exponen- 
tial increase of the scale factor. We can establish a useful general criterion for 
this asymptotic behaviour to occur for general a(t). Suppose that as t — > oo 
both sides of equation (H) tend to a constant (which may be equal to zero). 
Thus 



(^a 3 )' = A 

. ftdt f dt „ 

as t — > oo if 



A I *4 - (46) 



B 



Then for consistency we also require, as t — > oo, that 

ATexp[-2C] -> A 

so A cannot be chosen to be zero. Thus in all cosmological models for which 



(HH) and ([47]) hold there will be a self-consistent asymptotic solution as t 



oo of the form 



Vox P [-2(1 I** + b[% + C 



a° a 



where B,C are positive constants. We see that the n > 2/3 and de Sitter 
cases satisfy the conditions ( pEUD and ( fiTf ) hence ip is asymptotically constant. 
The dust (n = 2/3) case, for which tp — > oo as t —>■ oo satisfies (^) but 
violates 
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5 Conclusions 



Using a phase plane analysis we have studied the cosmological evolution of a 
time-varying fine-structure 'constant' a(t) = exp[2ip], in the BSBM theory. 
We have considered the cases created by power-law evolution of the expansion 
scale factor of the universe. We have shown that in general a increases with 
time or asymptotes to a constant value at late times. We have found a 
new exact solution for the case of a universe dominated by a stiff fluid or 
massless scalar field. We have also found general asymptotic solutions for all 
the different and possible behaviours via the analysis of the critical points of 
the system that determines the evolution of a. In particular, we have found 
asymptotic solutions for the dust, radiation and curvature-dominated FRW 
universe which also generalise the asymptotes found in 0. These solutions 
correspond to late-time attractors that describes the ip and a evolution in 
time and will enable a more detailed analysis to be made of the fit between 
theoretical expectations of varying a theories and observations of relativistic 
fine structure in atoms at high redshift. 
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Appendix: General analysis of the phase plane bifurcations 

In previous sections we have analysed the ip evolution equation (|8|) for a 
range of variables which are physically realistic and correspond to expanding 
universes. We will now analyse the whole range for variables of the system 
(|37|). As before we see there are two critical points in the (w,v) plane, at 
(0,-1 + f ) and ((1 - 3n)(|* - 1),0). Linearising fl37|) about (w Cl ,v Cl ) = 
(0,-1 + tP) and (w C2 ,v C2 ) = ((1 - 3n)(^ - 1),0) we obtain the following 
characteristics matrices: 
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Mi 



|orn= §. 



2-3n \ /O (3n-l)(2-3n) 

1 1 - 3n J 2 " V 1 1 - 3n 

The characteristic matrices are non singular except when n = 
In the non-singular cases the critical points will be simple and the system 
defined by these differential equations is structurally stable JJT], and there 
will be no 'strange' chaotic behaviour outside the neighbourhood of the crit- 
ical points. Hence, the linearised system will have the same phase portrait 
as non-linearised one in the neighbourhood of the critical points. 

The evolution, with respect to changing n, of the signs of the determinant 
and the trace of these two matrices is given in the table. This show us that 
there are always two critical points in our system, an unstable saddle and an 
attractor (which changes from a spiral to a node). 



n 


Critical Points (w c , v c ) 


((3n-l)(l-f),0) 


(0,1- f) 


(-oo; ±) 


Saddle Point (non-physical) 
det Mi < 0, Tr Mi > 


Unstable Node 
det M 2 > 0, Tr M 2 > 


l 

3 


Origin 
det Mi = 0, Tr Mi = 


Axis 

det M 2 = 0, Tr M 2 > 


(-■-) 

\3' 2/ 


Stable Spiral 
det Mi > 0, Tr Mi < 


Saddle Point 
det M 2 < 0, Tr M 2 > 


1 

2 


Stable Spiral 
det Mi > 0, Tr Mi < 


Saddle Point 
det M 2 < 0, Tr M 2 = 




Stable Spiral 
det Mi > 0, Tr Mi < 


Saddle Point 
det M 2 < 0, Tr M 2 < 


3 

5 


Stable Spiral (node) 
det Mi > 0, Tr Mi < 


Saddle Point 
det M 2 < 0, Tr M 2 < 


(hi) 
^5> 3^ 


Stable Node 
det Mi > 0, Tr Mi < 


Saddle Point 
det M 2 < 0, Tr M 2 < 


2 

3 


Stable Axis 
det Mi = 0, Tr Mi < 


Stable Axis 
det M 2 = 0, Tr M 2 < 


(|;oo) 


Saddle Point (non-physical) 
det Mi < 0, Tr Mi < 


Stable Node 
det M 2 > 0, Tr M 2 < 



The cases n 



or n 



3' 



where the determinant of the characteristic 
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matrixes vanishes, lead to a bifurcation of codimension 1, in particular, of 



Saddle-Node type pOfl , since they correspond to points where the determi- 
nants of the characteristic matrices change sign, detMi = detik/2 = 0. At 
these values of n the nature of the system will change. Cosmologically, these 
points represent a change in the behaviour of the time evolution of the fine 
structure 'constant' as can be seen from the figures: |], f|, |l], 0, ||, |3], which 
display the time evolution of ip. When n starts to grows from —00 to \ the 
two critical points slowly converge at n = | . For example, in the n = case 
where we may without loss of generality set a = 1 , equation (|8]) has the exact 
solution 

a = exp[2^] = A~ 2 cosh 2 [AN 1/2 (t + t )] 

where A, to are constants. This is an unrealistically rapid growth asymptoti- 
cally , t/> oc t, caused by the absence of the inhibiting effect of the cosmological 
expansion. The case of n = 1/4 is shown in Figure 6, which shows the phase 
space trajectories and the evolution of ip vs. Int. 




Figure 6: Numerical plots of the phase space in the (w,v) coordinates, and the tp 
evolution with x = log{t) for n = j. The '+ ' sign is a saddle point and the square 
is an unstable source. 

At n = I we are in the situation where the two critical points collapse 
into a unique one at the origin, creating a saddle-node bifurcation and a 
concomitant change in the behaviour and evolution of ip. As n keeps growing 
the single critical point splits into two critical points again. They move apart 
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until the radiation value is reached, n = \ (Tr M 2 = 0). In this case, ip 
is a asymptotically monotonic growing function of time, with some small 
oscillations near the Planck epoch. However, note that in our universe the 
asymptote giving an increase of ip behaviour with time is never reached before 



the dust-dominated evolution takes over 10 



As the universe evolves to the dust-dominated epoch, and n approaches 
the intermediate behaviour n = |, the two critical points start to coalesce 
again into a single point. When n = § is reached, ip becomes a strictly mono- 
tonically growing function of time. When n reaches the value corresponding 
to a dust-dominated universe, n — |, another saddle-node bifurcation occurs. 
The two critical points collapse into a single one. Again there will be change 
in the behaviour of ip for larger values of n. Accordingly, when n > 2/3, 
the two critical points reappear once again and ip becomes asymptotically 
constant in value. 

Notice, that although a bifurcation is something that 'spoils' the smooth 
behaviour of a system, in our case, that won't happen, due to the physical 
constraints of our variables. In reality due to those constraints, the physical 
system will never 'feel' the abrupt change at n = | and n = |. This is also 
due to the fact that the attracting critical point always lies in the physical 
range of the variables, while the unstable one disappears form the physical 
system when the bifurcations occur, as can be seen from the phase plane 
plots. 
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